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ABSTRACT 


The running skyline is a popular cable logging system being 
used in many parts of North America today. Proper application of 
this system requires investigation of its adaptability to the terrain 
being logged. This paper presents a digital computer program to 
determine running skyline load-carrying capabilities to aid the 
logging layout designer with this planning task. The Fortran IV 
program, the details of input, and the interpretation of output are 
discussed. 
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NOMENCLATURE FOR THE RUNNING SKYLINE COMPUTER PROGRAM.... LiL 
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INTRODUCTION 


The running skyline is a cable sys- 
tem widely used for logging in the timber 
harvesting industry. It is a system of 
two or more suspended, moving lines, 
generally referred to as a main and haul- 
back, that, when properly tensioned, will 
provide lift andtravel to the carriage. The 
basic system is illustrated in figure 1. 


A logging system designer must 
know what payload a logging system can 
carry over a givengroundprofile. Recog- 
nizing this requirement, Mann presented 
a procedure for determining the load- 
carrying capability of the running skyline 
and discussed the mechanics of some of 
the system's configurations. Mann's 


1/ Charles N. Mann. Mechanics of running sky- 
lines. USDA For. Serv. Res. Pap. PNW-75, 1lp., 
illus., 1969. Pac. Northwest For. & Range Exp. Stn., 
Portland, Oreg. 
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procedure was an extension of the 
graphical-tabular method discussed in 
the "Skyline Tension and Deflection 
Handbook, " 2/ 


Mann's procedure provides a 
straightforward approach to determina- 
tion of payload capability of the running 
skyline at midspan. Unfortunately, the 
method can be time-consuming, generally 
proceeding at a rate of 10 to 20 payload 
determinations per man-day. This rate 
may not be acceptable when a large num- 
ber of skyline roads are to be designed. 
A more efficient approach employs a 
desk-top computer/plotter. 3/ However, 


2/ Hilton H. Lysons and Charles N. Mann. Sky- 
line tension and deflection handbook. USDA For. Serv. 
Res. Pap. PNW-39, 41 p., illus., 1967. Pac. North- 
west For. & Range-Exp. Stn., Portland, Oreg. 


3/Ward W. Carson, Donald D. Studier, and 
Hilton H. Lysons. Running skyline design with desk- 
top computer/plotter. USDA For. Serv. Res. Note 
PNW-153, 21p., illus., 1971. Pac. Northwest For. 
& Range Exp. Stn., Portland, Oreg. 
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Figure 1.--Running skyline. 


when the layouts involve many skyline roads 
and desk-top facilities are not available, 
it is expeditious to employ a high-speed 
digital computer for direct computation 
of payload capability. A computer pro- 
gram designed for this purpose is pre- 
sented and discussed in this paper. 


This computer program is an ana- 
lytical tool that canbe used by the timber 
resource planner for determining the 
feasibility of the running skyline system. 
It has been prepared to provide the logging 
systems designer with a tool for comput- 
ing the vertical load-carrying capability 
of a grapple-rigged running skyline. The 
program requires information about the 
ground profile over which the system is 
operating, the equipment used, and the 
tension to be maintained in the haulback 
line during the yarding process. The 
method of solution is essentially that 
discussed by Carson and Mann.4/ The 
program was written in standard Fortran 
IV language for the CDC 6400 machine. 
It has also been prepared to operate on 
the CDC 3100. An attempt has been made 
to keep the language standard enough that 
it will operate on most machines having 
a Fortran IV compiler. 


PROGRAM LIMITATIONS 


The designer using this program 
should be aware of some of the limita- 
tions of the results. The program has 
been designed to compute the load- 
carrying capability of a grapple-rigged 
running skyline system. These computa- 
tions assume that an operating tension 
exists in the haulback line, that the load 
is suspended vertically below the carriage, 
and that the other force required to hold 


4/ Ward W. Carson and Charles N. Mann. A 


technique for the solution of skyline catenary equations. 


USDA For. Serv. Res. Pap. PNW-110, 18 p., illus., 
1970. Pac. Northwest For. & Range Exp. Stn., Port- 
land, Oreg. 


the log and carriage in position is provided 
by the main line only. The main line 
tension is limited to a maximum value; 
therefore, either the haulback operating 
tension or the main line maximum tension 
may decide the load-carrying capability. 


In an actual situation, it is quite 
possible that one end of the log will be 
dragging on the ground during the yard- 
ing process. In this case, the full weight 
of the log will not be felt by the running 
skyline system. This case differs in two 
significant ways from the situation where 
the log is suspended clear of the ground: 
and hangs directly below the carriage: 
the vertical force on the carriage is less 
than the log weight, and there is a hori- 
zontal force on the carriage due to the 
log dragging. This case is not treated 
by the program discussed in this paper. 
For such situations, this program can 
only provide the designer with some quan- 
titative feel for the system's capability. 


Caution must be exercised in apply- 
ing the results of this program to a running 
skyline system whichuses a slack-pulling 
carriage. When a Slack-pulling carriage 
is used, it is possible to design the equip- 
ment so that the yarding force is shared 
by both the main line andthe slack-pulling 
line. Since this program assumes that 
only the main line bears this force, it 
does not apply in a shared situation. In 
those cases where the main line does bear 
the force and the log is hanging below the 
carriage, the results are applicable. 


Therefore, caution should be exer- 
cised when the results of this program are 
used for any running skyline configuration 
other than the grapple carriage supporting 
a completely suspended load, or a slack- 
pulling carriage supporting a suspended 
load with negligible tension in the slack- 
pulling line. When other situations exist, 
this program is suggested merely as a 


tool to give the designer a conservative 
estimate for the capabilities of the running 
skyline system. 


PROGRAM DESCRIPTION 


The basic data needed as input for 
the grapple-rigged running skyline pro- 
gram consist of equipment specifications 
and individual skyline road specifications. 
More than one set of skyline road speci- 
fications may be used with each set of 
equipment specifications. 


The equipment specifications include 
the weight of the carriage; the carriage 
clearance above ground; the headspar and 
tailspar heights; and the diameter, weight 
per foot, breaking strength, and safety 
factor of the haulback and main line cables. 
The weight and breaking strength for a 
given cable diameter can be obtained 
from "Skyline Tension and Deflection 
Handbook" (see footnote 2). 


The individual road specifications 
are made by describing the headspar and 
tailspar locations, the inner and outer 
yarding limits, and the range and eleva- 
tion of enough profile points to character- 
ize the ground under the skyline. The 
following convention was adopted for these 
data (fig. 1): 

1. The yarder is always on the left 

at spar 1. 

2. The terrain points are expressed 

as coordinates: the abscissa x 
denotes the horizontal distance 
measured from spar 1, or from 
some reference point to the left 
of spar 1, in feet; the ordinate y 
denotes the elevation of the point, 
in feet. 

All of the above specifications are read 

in, and values are computed for internal 

variables in sections 1 and 2 of the main 

program. 


In section 3, the geometric parame- 
ters, such as slope (S), span (L), the 
difference in elevation of the top of the 
spars (H), and the carriage's vertical 
(DY1) and horizontal (D1) position are 
determined from the input data. 


In section 4 of the main program, 
each terrain point is checked to deter- 
mine if it intersects the chord of the 
skyline and if it is within the yarding 
limits. These checks are performed 
before any skyline computations are 
made. [If a terrain point is found to inter- 
sect the skyline chord, the computations 
for that skyline road are terminated, and 
a diagnostic statement is written. Sky- 
line computations are performed only on 
terrain points that are within the yarding 
limit. The yarding limits can be located 
at the headspar, tailspar, or anywhere 
within the skyline span. 


The line tensions and payload are 
computed in section 5 of the main pro- 
gram. Before these quantities can be 
determined, the geometric configuration 
of the skyline must be known because the 
analysis depends on the relative positions 
of spar 1, spar 2, and the carriage. The 
allowable tensions of the haulback line 
and the main line determine in which line 
the tension is critical, and the critical 
tension always occurs at the highest 
point reached by a given line. 


If the yarder is at the upper end, 
as shown in figure 2, the carriage will 
always be below the top of spar 1, and 
the critical tension will occur at spar 1 
(Type = 6), either in the main line or in 
the haulback line (Type = 8) depending on 
their relative sizes. If the yarder is at 
the lower end, one of the three possibili- 
ties depicted in figures 3, 4, or 5 can 
occur. The critical tension may occur 
in the haulback line at spar 2 (Type = 7), 
regardless of the position of the carriage, 


Critical tension (haulback line) 
Critical tension (main line) 


ak Haulback line 
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Figure 2.--Yarder at upper end. Critical tension occurs at spar l, 
either in the haulback line (Type = 6), or in the main line 
(Type = 8). 
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Figure 3.--Yarder at lower end. Critical tension occurs at spar 2 
in the haulback line (Type = 7). 
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Figure 4.--Yarder at lower end. Carriage below top of spar l. 
Critical tension occurs at spar 1 in the main line (Type = 8). 
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Figure 5.--Yarder at lower end. Carriage above top of spar l. 
Critical tension occurs at the carriage in the main line 


(Type = 9). 


as shown in figure 3. However, if the 
line sizes are such that the critical ten- 
sion occurs in the main line, then the 
position of the carriage will determine 
its location. If the carriage is below the 
top of spar 1, the critical tension will 
occur at spar 1 (Type = 8), as shown in 
figure 4. _If the carriage is positioned 
above the top of spar 1, the critical ten- 
sion will occur at the carriage (Type = 9), 
as shown in figure 5. 
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After all the geometric parameters 
of the skyline road and the type are deter- 
mined, control of the program is trans- 
ferred to subroutine CONVG. This sub- 
routine controls the iterations for solving 
the catenary equations (see footnote 4). 
The flow diagram illustrating the iterative 
procedure is shown in figure 6. The 
initial estimates for the catenary parame- 
ters m,, m,, and m3 are found by solv- 
ing for horizontal tensions in the force 
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Figure 6.--Macro flow diagram: 


Main program. 


balance equation. This is done in sub- 
routine FB3. These initial estimates for 
the m's are then substituted in the cate- 
nary equations in subroutine CATGEO 
and tested for convergence in subroutine 
TEST3. The convergence equations are 
derived from the boundary conditions 
related to the type of skyline configura- 
tion. If the convergence functions are 
found to be within the tolerance set in 
TEST3, the control of the program is 
returned to the main program with the 
solution. If the convergence functions 
do not fall within the tolerance, new line 
lengths and moment arms are calculated 
in subroutine CATARM for use in FB3 

to provide new estimates for the m's, 
The process is then repeated until a solu- 
tion is found. 


At this point, all conditions have 
been satisfied and the program execution 
could be terminated. However, one more 
step is taken. Determination of the 
unloaded carriage clearance at the given 
terrain point--that is, the distance from 
the ground to the carriage when the pay- 
load is equal to zero--is done in sub- 
routine DEFL3. Determination of the 
unloaded carriage clearance involves Ay 
as an unknown. The iterative procedure 
of this routine employes the secant meth- 
od. This method requires two initial 
estimates of the unknown Ay, plus a 
third estimate that is calculated from 
the secant formula. The first estimate 
of Ay is taken from the solution of the 
skyline problem along with the calculated 
payload. The second estimate of Ay is 
taken as a percentage of the first, and a 
new payload is calculated by entering 
subroutine CONVG. If the payload is 
found to be unequal to zero, the secant 
iteration loop is entered, a new Ay is 
determined from the secant formula, and 
subroutine CONVG is again entered to 
calculate a new payload. If this new pay- 
load is found to be unequal to zero, the 


secant formula is again used to calculate 
anew Ay, and the process is continued 
until the payload is found to equal zero. 


INPUT 


Input consists of the title, the equip- 
ment data, and the geometric data. Sample 
data forms are shown in figure 7. The 
title, which is contained on two cards, is 
read in an alphanumeric format. The 
first card contains one field of 70 columns 
for the project title and one field of 10 
columns for the date. The second title 
card contains four fields: 10 columns 
for the region, 20 columns each for the 
Forest and District, and 30 columns for 
the location of the project. 


The equipment data are also con- 
tained on two cards. The first card 
consists of eight fields of 10 columns, 
read by the Fortran 8F10. 1 format. 
Decimal points must be punched. The 
cable diameter, weight, and breaking 
strength for both the haulback line and 
main line, plus the safety factor and 
carriage weight, are contained on this 
card. The second equipment data card 
specifies the number of skyline roads 
that will accompany the set of title and 
equipment data, and the input type, read 
by the Fortran 25, 155 format. 


The input type is coded zero if the 
geometric data are furnished by the de- 
signer in the form of x and y coordinates. 
If the geometric data are to be taken 
from aerial photographs, input type is 
coded 1. The procedure for submitting 
aerial photos is discussed in Research 
Note PNW-132.3/ 


5/ Ward W. Carson, Donald D. Studier, and 
William M. Thomas. Digitizing topographic data for 
skyline design programs. USDA For. Serv. Res. Note 
PNW-132, 13 p., illus., 1970. Pac. Northwest For. 
& Range Exp. Stn., Portland, Oreg. 


RUNNING SKYLINE 


Forest: District: Location: 


Equipment Data 
ie a eS HAULBACK Ss ines (00 a |e ase MATING N ae) 
Diameter Weight Strength Diameter Weight Strengt Factor Weight 
XX .XX XXX XX XXXXXX XX XX XX XX .XX XXXXXX XX XXX XXXX «XX 


Format of Geometric 80 
Input Data 


Number of 
Skyline Roads 


Spar 1 Haulback line: 


Main line: Diameter (inches) 
Diameter (inches) Weight per foot (pounds) 
Weight per foot (pounds) Breaking strength (pounds) 
Breaking strength (pounds) 
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Weight (pounds) 
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x = Horizontal distance from left 
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Clearance = Distance from the car- 
riage to the ground (feet). 
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Figure 7.--Input data form. 


One set of geometric data describes 
an individual skyline road. The first geo- 
metric data card contains the skyline road 
number and is read by the Fortran 125 


set of geometric data. When all sets of 
geometric data are read, the program 
will proceed to the next set of title and 
equipment data with its related geometric 


format. The second geometric card con- data. 
sists of nine fields of five columns each, 
read by the Fortran F5.0 format. This 

OUTPUT 


card contains the carriage clearance, 
height and location of the left and right 
spars, and the inner and outer yarding 
limits. The terrain data cards contain 
16 fields of five columns each, read by 
the Fortran 16F5.0 format, and the values 
of distance (x) and elevation (y). There 
can be more than one terrain data card; 
the only requirement is that the last card 
be left blank. This blank card marks 

the end of the data for that skyline road, 
and the program will proceed to the next 


The output is illustrated in figure 8. 
The title and equipment data are printed 
in an orderly manner for easy reference. 
This is followed by the data for each sky- 
line road, which consist of the input data 
and the solution to the problem. 


The first section of the output con- 
sists of the title and equipment data that 
were furnished by the designer. It includes 
the project title, date, and location; the 


LEWIS RIVER SKYLINE SALE 4/10/70 
REGION 6/ FOREST  GIFFORD PINCHOT/ DISTRICT LEWIS RIVER/ LOCATION LEWIS RIVER 
CABLE DIAMETER LINE WEIGHT BREAKING STRENGTH 
HAULBACK 75 1.04 58800 .00 
MAIN LINE 1.00 1.85 103400 .00 
SAFETY FACTOR = 3.00 
CARRIAGE WEIGHT = 5000.00 KK KIKI KK KKK KKK KKKK KE RKEKE 
ROAD NUMBER = 1 LOAD CLEARANCE = _—20.0 
SPAR HEIGHT STATION ELEVATION 
HEAD SPAR 50.0 0.0 4826 .0 
TAIL SPAR 20.0 1500.0 4734.0 
YARDING LIMITS 0 1500 
SKYLINE SPAN = 1500.0 SLOPE = -8.1 PERCENT 
RUNNING SKYLINE OUTPUT AT ALL TERRAIN POINTS 
“TOT. te a es a CABEE TIEN SI ONS saan UNLOADED 
PAYLOAD HAULBACK HAULBACK MAIN LINE CARRIAGE 
BUST NG ELEVATION LBS. DEFLECTION YARDER TAILSPAR YARDER CLEARANCE 
179.0 4792.0 5563.9 3.30 19600.0 19473.2 21852.7 42.5 
332.0 4771.0 2026.8 3.87 19600.0 19473.1 20728 .4 33.5 
510.0 4740.0 1836.0 4.97 19600.1 19473.1 20507 .8 35.9 
701.0 4702.0 3230.1 6.47 19600.1 19473.1 20436 .6 51.2 
891.0 4637 .0 8574.0 9.77 19600.3 19473.2 20238 .0 102.1 
996.0 4570.0 15541.8 1967 19600.4 19473.3 19203.2 165.0 


Figure 8.--Output format. 


diameter, weight, and breaking strength 
of the haulback and main lines; the safety 
factor; and the weight of the carriage. 
These data may be common to more than 
one skyline road. 


The second section of the output 
applies to one particular road. The road 
number, load clearance, head spar and 
tail spar heights, station, elevation, and 
yarding limits are data that were fur- 
nished by the designer. The program 
calculates the skyline span, which is the 
horizontal distance between the head spar 
and tail spar, and the percent slope of 
the skyline chord. 


The running skyline capability is 
calculated for each terrain point which 
falls within the yarding limits. [If the 
terrain point is outside the yarding limit, 
a diagnostic statement is printed. 


A discussion of the columns listed 
under "Running Skyline Output at all 


Terrain Points" (fig. 8) follows: 


The first two columns, horizontal 
distance and elevation, describe the 
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terrain point and are furnished by the 
designer or taken directly from the aerial 
photographs (see footnote 5). 


The calculated payloads, listed in 
the third column, are the payleads that the 
skyline will support when the carriage is 
positioned at the required clearance 
above the terrain points, 


Deflection is the vertical distance 
from the chord to the skyline, expressed 
as a percent of the span length. This 
program calculates the deflection in 


_percent, at the terrain point in question, 


when the skyline is supporting the calcu- 
lated payload. 


The next three columns list the 
cable tensions that occur at the yarder 
and tailspar in both the main and 
haulback lines when they are supporting 
the calculated payload. 


Inthe lastcolumn, the unloaded car- 
riage clearance is listed for eachterrain 
point. This is the distance between the 
carriage and the ground when the carriage 
is unloaded and the working tension is 
maintained in the lines. 
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Intermediate quantity used in the computation of horizontal tension, 
pounds squared. 


Intermediate quantity used in the computation of horizontal tension, 
pounds squared. 


Coordinates of terrain points, determined from aerial photographs. 


Coordinates of locations of the head- and tailspars and the inner 
and outer yarding limits, determined from aerial photographs. 


Required carriage clearance, feet. 
Distance between skyline chord and specified carriage location, feet. 


Elevation of carriage when it is at the specified clearance above 
the terrain point, feet. 


Cable diameter of haulback line, inches. 

Elevation of skyline chord, feet. 

Input title storage. 

Deflection at given terrain point. 

Horizontal location of terrain point from spar 1, feet. 
Deflection, percent. 


Intermediate quantity used in the computation of horizontal tension, 
pounds squared. 


Intermediate quantity used in the computation of horizontal tension, 
pounds squared. 


Vertical distance from spar 1 to carriage when it is located at the 
specified clearance distance above the terrain point, feet. 


Moment arm of the weight of line segment 1 from the carriage, feet. 
Moment arm of the weight of line segment 2 from the carriage, feet. 


Moment arm of the weight of the main line from the carriage, feet. 
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Elevation of spar 1, feet. 

Elevation of spar 2, feet. 

Convergence functions for the various ITYPEs. 

Difference in elevation between top of spar 1 and top of spar 2, feet. 


Horizontal component of tension in line segment 1 of the haulback 
line, pounds. 


Horizontal component of tension in line segment 2 of the haulback 
line, pounds. 


Horizontal component of tension in the main line, pounds. 

Height of spar 1, feet. 

Height of spar 2, feet. 

Indicator that functions have converged. 

Indicator, if equal to 1, causes intermediate output to be printed. 
Indicator for ITYPE. 

Indicator noting the first terrain point inside the inner yarding limit. 
Indicator noting the last terrain point inside the outer yarding limit. 
Indicator used in routine accounting for yarding limits. 

Indicator used to eliminate some calculations after the first iteration. 


Indicator noting if the input data is two-dimensional or three- 
dimensional. 


Horizontal distance between spar 1 and spar 2, feet. 
Horizontal location of spar 2, feet. 

Ratio of DY and D. 

Ratio of El and D. 

Ratio of ES and DS, 


Ratio of (DY-H) and (L-D). 


ST5 


ST6 


TA 


TAS 


TBSKY 


TBSNUB 


TITLE1 


TIT LE2 


TOL 


v1 


V2 


VS 


WC 


WL 


Ratio of E2 and (L-D). 

Ratio of SYS and DS. 

Allowable tension in haulback, pounds. 

Allowable tension in main line, pounds. 

Breaking strength of haulback line, pounds. 

Breaking strength of main line, pounds. 

Input title storage. 

Input title storage. 

Tolerance allowed in the convergence function. 

Vertical component of tension in line segment 1 of the skyline, pounds. 
Vertical component of tension in line segment 2 of the skyline, pounds. 
Vertical component of tension in the main line, pounds. 

Weight of the carriage, pounds. 

Payload or weight of logs, pounds. 

Weight of segment 1 of haulback, pounds/foot. 


Weight of segment 2 of haulback, pounds/foot. 
Weight of main line, pounds/foot. 


Horizontal location of terrain point, feet. 

Geometric quantity used in the solution of catenary equations. 
Geometric quantity used in the solution of catenary equations. 
Geometric quantity used in the solution of catenary equations. 
Geometric quantity used in the solution of catenary equations. 
Geometric quantity used in the solution of catenary equations. 
Geometric quantity used in the solution of catenary equations. 


Elevation of terrain point, feet. 


Horizontal location of inner yarding limit, feet. 


Horizontal location of outer yarding limit, feet. 
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RUNNING SKYLINE PROGRAM 


C PROGRAM RSKY 


C ITYPE=6 HeGE.O 

ITYPE=7 HeLTeO 

ITYPE=8 DY-GEe0 

ITYPE=9 DYeLT.O 
MACHINE DATA 
TBSKY =BREAKING STRENGTH OF HAULBACK»s POUNDS 
TBSNUB=BREAKING STRENGTH OF MAIN LINEs POUNDS 
CDSKY =DIAMETER OF HAULBACK»s INCHES 
CDSNUB=DIAMETER OF MAIN LINE? INCHES 


WO1  =WEIGHT OF HAULBACK»s POUNDS/FOOT 
WOS -=WEIGHT OF MAIN LINEr POUNDS/FOOT 
SF =DESIGN SAFETY FACTOR 

WC =WEIGHT OF CARRIAGE? POUNDS 

HS1  =HEIGHT OF SPAR 1+ FEET 

HS2. =HEIGHT OF SPAR 2+FEET 


MISCELLANEOUS DATA 

YLAi=INNER YARDING LIMIT? FEET 
YL2=OUTER YARDING LIMIT» FEET 
ES1=ELEVATION OF SPAR 1+FEET 
ES2=ELEVATION OF SPAR 2e°FEET 
SS1=STATION OF SPAR ir FEET 
SS2=STATION OF SPAR 2° FEET 
IRN=SKYLINE ROAD NUMBER 


Cy GEOMETRY 
S =SEOPE 
L =SPAN 
DI =STATION 
DYI=DISPLACEMENT FROM TOP OF SPAR i TO CARRIAGE 
H =DISPLACEMENT FROM TOP OF SPAR 1 TO TOP OF SPAR 2 
Y =TERRAIN POINT ELEVATION 
X =TERRAIN POINT STATION 
C1 =REQUIRED CARRIAGE CLEARANCE 
DIMENSION X(50)»Y(50) + DI(50) »DYI(50) » CACL( 50) » WP1(50) »DF (50) «CL (50 


ALL TN) FREER 


NANANANANANAN DONANIDIANIANAANANIADAIAAANANDNDIANANANAAANNADA 


1) 

DIMENSION AX(50) #AY(50) # AZ( 50) » TA( 50) 

DIMENSION BX(50) eBY(50) + BZ(50) + IB( 50) 

DIMENSION TITLE1(16) » TITLE2(16) e DATE( 2) 

DIMENSION T1(50) eT2(50) » TS( 50) 

REAL LeLireL2*+LSeMO1*»MO2* MOS 

COMMON IPRINTe* ITER 

COMMON TArTASeWO1’ WO2*+WOSeL»e Se DEFL?s De DY DSt He DYSe WCe Cle HS1#% HS2 


COMMON 
COMMON 


LieL2eLSeE1eE2°ES 
MO1» MO2*MOS+* XAMO1* X1MO1» XAMO2¢ X2MO22 XAMOS? X1MOS 


COMMON wWLeSG 
INTEGER OT 
C SECTION 1 *¥**** INITIALIZE 
IN=60 
OT=61 
NP=62 


14 


IPRINT=i1 
IPRINT=0 
(c 
C SECTION 2 ****xX INPUT = OUTPUT 
1 CONT INUE 
READ (IN*e1005) (TITLEL( I)» I=1014) »DATE(1) »DATE(2) 
GO TO (5000r2) EOFCKF (IN) 
2 READ (INe1005) (TITLE2(1I) » I=1°16) 
/ ISET=0 
| WRITE (OT*#2001) (TITLE1( I) » I=1912) »DATE(1) »DATE(2) 
WRITE (OTe2002) (TITLE2( I)» I=1°16) 


(c 
C . EQUIPMENT SPECIFICATIONS 

READ (INe1000) CDSKY»wO1e TBSKY» CDSNUB?e WOSe TBSNUB»s SFe WC 
SF=5.0 
WRITE (OT»2004) CDOSKY»WO1» TBSKY» CDSNUB?p WOSs TBSNUBp SFe WC 
TA=TBSKY/SF 
TAS=TBSNUB/SF 
woO2=wOl 
READ(IN»e 1001) NOR»eSCL»e INTYP 


(E 
C ‘SKYLINE ROAD SPECIFICATIONS 
3 IF (INTYPeEG.0) GO TO 9 
(C 
C AUTO=TROL DATA 
IF(SCL.«EQe020) SCL=1.0 
READ (IN»1004) IRNeC1»HS1»HS2e TECARDe (BX( 1) »BY( 1)» BZ(1I)»+ IBC I) » 
1I=194) 
Ni=1 
uy N2=N1i+3 
READ (IN*#1006) (AX( I) e AY( I)» AZ( I)» TACIT)» T=Nip N2) 
NL=Ni+4 
IF (IA(N2).EQe2) GO TO 4 | 
5 IF (IA(N2) -EQ.3) GO TO 6 
N2=N2=1 
GO TO 5 
6 N=N2 
SS1=SQRT( (BX (1) =AX(1) ) **2+4+(BY(1) =AY(1) )**2) *SCL 
ES1=BZ(1) 
YL1=SQRT( (BX(2)=AX(1) ) ®*24( BY (2) AY( 1) ) **2) *SCL 
YL2=SQRT( (BX(3)~=AX(1) ) ®*24+(BY (3) =AY(1) ) **2) *SCL 
SS2=SQRT ( (BX (4) AX (1) ) **2+ (BY (4) AY (1) ) *¥*2) *SCL 


ES2=BZ (4) 

WRITE (OT*#2005) IRNe C1» HS1*°SS1irESi*HS29SS2eES2eYL1ir YL2 
DO 7 I=19N 

X( I) =SQRT( (AX (T)=AX(1) ) ®*24+ (AY (1) AY (1) ) #2) *SCL 

Y( I) =A2(T) 


7 CONTINUE 
WRITE(NP*®1001) IRN 
WRITE (NP»1003) Ci»HS1i*SS1*ES1*HS29SS2eES2e YL1» YL2 
WRITE(NP? 1002) (X(T) e YC 1)» T=19N) 
WRITE( NP» 7000) 
GO TO 30 

C 

C TERRAIN DATA=NOTE* ONE SET OF BLANKS REQUIRED AFTER LAST SET OF TERRAIN 


9 READ (IN»e1001) IRN 
GO TO (1°91) EOFCKF(IN) 
91 READ (IN®1003) Ci»HS1*»SST*+ES1i» HS2*°SS2°ES2e YL1ir YL2 
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WRITE (OT®2005) IRNeC1*eHSirSSir ES1*®HS2*SS2eES2e YL1ir YL2 
N1=1 
N2=8 
10 READ (IN*®1002) (X( 1)» Y(1)» T=N19N2) 
DO 11 I=N1°N2 
IF (Y(I)-EQe-0-0) GO TO 12 
11 N=I 
N1=N1+8 
N2=N2+8 
GO TO 10 
12 IF(NeNEe(Ni-1)) READ (IN? 1001) 
Cc 
C SECTION 3 ****k* ESTABLISH GEOMETRY 
30 L=SS2-Ss1 
S=((CEST HSI i= (ES2HHS2))) 7 
SP=-100-0*S 
H=S*L 
WRITE (OT»2006) LeSP 
DO 32 I=19N 
DI(1I)=X(I)-SS1 
DYI( I) 2(ESI+HS1)-Y( 1) =C1 
32 CONT INUE 
C CHECK FOR YARDING LIMITS OUTSIDE OF THE ANCHOR POINTS 
IF(YL1eLTeSS1e¢OR.YL2eGT«SS2) WRITE( OTs 2019) 
IF(YL1-LTeSS1.eO0ReYL2eG6T.«SS2) GO TO 100 
C 
C SECTION 4 ***k* ESTABLISH TERRAIN POINTS 
C 
IL1=1 
IL2=N 
ISET=0 
C CHECK YARDING LIMITS 
DO 40 I=19N 
TE OX) eee Et) GO) TON as 
IF (X(I)eGEeYL2) GO TO 46 
IF (X(T) eLTeSS1-¢OReX(1)eGTeSS2) GO TO 421 
GO TO 47 
45 CONTINUE 
IL1=I+1 
WRITE (OT*#2010) X(I)*Y(T) 
GO TO 40 
46 CONT INUE 
ILT=I-1 
TL2=MINO(IL2e ILT) 
WRITE (OT*#2010) X(T) Y(T) 
GO TO 40 
C 
47 CONT INUE 
C CHECK TERRAIN FOR CHORD INTERSECTION 
CHEL=ES1+HS1=DI(I)*S 
CAEL=Y(1I)+C1 
CACL( 1) =CHEL-CAEL 
DF( IT) =CACL(I)/L*100.0 
DEFL=CACL(I)/L 
IF (CACL(I).LE.0.0) ISET=1 
GO TO 40 
421 WRITE(OT®2017) X(T) Y(TI) 
GO TO 100 
40 CONTINUE 
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C 
GS 
C 


50 
31 


52 


53 
S54 


60 
61 


99 


100 


5000 
1000 
1001 
1002 
1003 
1004 
1005 
1006 
2001 
2002 


LE GrSE lke EQ 1s)" GO" 10) 61 
ECTION 5 **** RUNNING SKYLINE COMPUTATIONS 


DO 60 I=IL1rIL2 
D=DI(T) 
DS=D 
DY=DYI(T) 
DYS=DY 
DETERMINE ITYPE 
IF (HeGEe-0-0) GO TO 50 
ITYPE=7 
GOnTO 52 
ITYPE=6 
CALL CONVG(ITYPE) 


CHECK DY 
IF(DYeGEe0-0) GO TO 52 
TDIFF=TAS=wOS*MOS* COSH( XAMOS) 
IF (TDIFFeGEe0-0) GO TO 54 
ITYPE=9 
GO TO 53 
TDIFF=TAS-wOS*MOS* COSH( X1MOS) 
IF (TDIFFeGE-0.0) GO TO 54 
ITYPE=8 
CALL CONVG(ITYPE) 
WP1(1I)=WL 
T1( 1) =[MO1*WwO1* COSH(X1M01) 
T2( 1) =MO2*WO2*COSH( X2M02) 
TS( I) =MOS*wWOS* COSH( X1MOS) 
CAEESDERESCITYPE) 
CL( 1) =HS1+ES1-Y ( 1)- (D*¥S+DEFL*L) 
CONT INUE 
CONT INUE 
IF (ISETeEQe1) WRITE (OT* 3000) (X(T) Y( I)» CACL( 1)» T=IL1» IL2) 
IF (ISET»EQe1) GO TO 100 
CONT INUE 
WRITE ( OT» 2007) 
WRITE (OT*#3001) (X( 1) Y(I) eo WP1(1) + DF( I) T1(T) TACT)» TS(T) CLC I)» I= 


Pee TE) 

CONT INUE 

GO TO 3 

RETURN 

FORMAT ( 8F 10.0) 

FORMAT (20X* 150 5X*eF1060° 35X09 15) 

FORMAT (16F5.0) 

FORMAT (16F5.0) 

FORMAT (I503F3-001194(3F5-0911)) 

FORMAT (16A5) 

FORMAT (15X04(3F5.0011) ) 

FORMAT (1H1»920X»912A59 20X» 2A5) 

FORMAT (8H REGION »2A509H/ FOREST #4A5e11H/ DISTRICT #4ASe 
1 12H/ LOCATION » 6A5/) 


2004 FORMAT (1H #20X015H CABLE DIAMETER? 5X012H LINE WEIGHT» 8X018H BREA 


LING STRENGTH/9H HAULBACKe 5X0 3F20¢2/9H MAINLINE? 5X0 3F20-2/15H S 
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2Y FACTOR=0F5e2/17H CARRIAGE WEIGHT=? F8-1//) 
2005 FORMAT (1X0 4/0 40X90 QOHKFEE KKH KEE KEK KEREEE 9 // 1 


3 13H ROAD NUMBER=» 15°20X016H LOAD CLEARANCE=9F5.1/23X012H § 


1PAR HEIGHT 10X+ 8H STATION» 9X¥*#10H ELEVATION/12H HEAD SPAR? 3F20.1/ 


212H TAIL SPAR*# 3F2061//3X9 1 4HYARDING LIMITS? 2F10-) 
2006 FORMAT (15HO SKYLINE SPAN=sF7.1*5X* 7H SLOPE=s F6e1*8H PERCENT//) 


2007 FORMAT (45H RUNNING SKYLINE OUTPUT AT ALL TERRAIN POINTS/» 
TAL 0 2a meme ee eee ee GCABEE TENSTONS ee » 7X» BHUNLOADED / 


2° 3XeLOHHORIZONT AL? 25X» THPAYLOAD? 26X# BHHAULBACK # 9X# BHHAULBACKs 9 Xp 


38HMAINLINE? 7X» BHCARRIAGE/# 4X» BHDISTANCEr 8X QHELEVAT ION? 11X*4HLBS,, 


46X* LOHDEFLECT IONe 12X* GHYARDER? 10Xe BHT AILSPAR? 10X* GHYARDER? 8X» 
59HCLEARANCE ) 

2010 FORMAT (2F9e1*#34H THIS POINT OUTSIDE YARDING LIMITS) 

2012 FORMAT (2F9e192(2F9.19F9.4) # F901) 

2013 FORMAT (2F9.19°2(2F9.19F904) »F9e1*#18H HAULBACK REQUIRED) 

2017 FORMAT(2F9e1*°75H THIS POINT OCCURS OUTSIDE OF THE ANCHOR POINTS, 


$SCHECK THE YARDING LIMITS. 
2019 FORMAT (/18X*61H THE YARDING LIMITS MUST NOT BE OUTSIDE OF THE ANC 


SHOR POINTS /) 
3000 FORMAT(79HOCARRIAGE WILL NOT CLEAR THE GROUND AT TERRAIN POINTS WI 


$TH NEGATIVE CLEARANCE 730HO STATION ELEVATION CLEARANCE // 
$(3F10.0)/) 


3001 FORMAT (F11¢192(5X0F1201) »5SXeFBe2" BXF120193( 5X0 F1201) /) 
7000 FORMAT(1H ) 


END 


SUBROUTINE CONVG(ITYPE) 


C ROUTINE TO CONTROL ITERATIONS BETWEEN FORCE BALANCE AND CATENARY 
C 


REAL LeLieLo2relLSeMO1*»MO2»MOS 

COMMON IPRINTe ITER 

COMMON TArTASe W012? WO2+WOSr Le Se DEFL? De DY? DSt He DYSe WCe C1 # HS1 2 HS2 
COMMON LieL2eLSreE1eE2°ES 


COMMON MO1*+MO2+MOSe XAMO1* X1MO1 * XAMO22 X2MO2¢% XAMOS ep X1MOS 
COMMON WLeSG 


C 
C START ITERATION LOOP 
C 
IND=1 
ITER=0 
10 ITER=ITER+1 
c 


C GET FORCE BALANCE ESTIMATES OF HORIZONTAL TENSIONS 
CALL FB3( IND? ITYPE) 
C 


C GENERATE CATENARY GEOMETRY 
CALL CATGEO 
C 


C CHECK CONVERGENCE 
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CALL TEST3S(ITYPE> I) 
IND=0 


IF (ITEReEQ210) RETURN 
een eleaOet) GO 10. 20 
RETURN 
20 CALL CATARM 
GO TO 10 ‘ 
END 


SUBROUTINE FB3( INDs ITYPE) 


REAL LeLieL2eLSeMO1» M02» MOS 
COMMON IPRINT?* ITER Y 
COMMON TAreTASeWO19*WO29eWOSe Le Se DEFL?e De DY e DSe He DYSe WCe C19 HS19 HS2 
COMMON LireL2eLSe Ele E20 ES 
COMMON MO1*MO2% MOS XAMO1 » X1MO1 » XAMO2 » X2MO02 0 XAMOSe XLMOS 
COMMON WLeSG 
INTEGER OT . 
OT=61 , 
C COMPUTE NECESSARY COEFFICIENTS y 
C INITIALIZE IF NECESSARY 
IF (INDeNEe1) GO TO 10 R 
C INITIALIZE LINE LENGTHS AND MOMENT ARMS 
L1=SQRT (D**2+DY**2) 
L2=SQRT( (DY=H) ** 2+ (L=D) **2) 
LS=SQRT ( DS**2+DYS** 2) C 
E1=D/2.0 
E2=(L=D) /2.0 
ES=DS/2.0 
10 ST1=DY/D 
ST2=E1/D Vlg 
ST3=ES/DS 
ST4=(DY=H) /(L=D) 
ST5=E27 (E=D) 
ST6=DYS/DS b 
C ESTABLISH LINE WEIGHTS 
R1=Li*wOl " 
R2=L2*w02 } 
RS=LS*WOS , 
TIDUM ="ITYPE = 5 
GO TO (15920°25*30)» IDUM 
C 
(c FOR HeGEe0 = COMPUTE PAYLOAD TYPE=6 
C 
ANS CONT INUE % 
C COMPUTE MOL 
ARG1=(R1I*ST1I*ST2) **¥2—=(12e04+ST1I**2) *( (RI*ST2) €*¥2-TA** 2) 4 
IF (ARG1.LT.0.0) GO TO 90 
H1=(=R1*ST1*ST2+SQRT(ARG1))/7(1-¢0+ST1**2) 
MO1=H1/WO1 ¥ 
c 
C COMPUTE M02 
V1I=R1I*ST2+H1*ST1 
DUM1L=(V1=R1) **2+H1**2 > 
ARG2=(R2*ST4*(ST5=1.0)) *¥2=(1¢04ST4H*K*2) *( (R2*(ST5~71.0) )**2=DUM1) 


IF (ARG2-LT-0.0) GO TO 90 | 
H2=(=R2*ST4*(ST5—1.0) +SQRT(ARG2))/7(1¢0+ST4H* +2) 
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C 
C 


C 
C 


C 
C 
C 


C 
C 


C 


C 


C 
C 
(SG 


25 CONT INUE 


C 


C 
€ 


20 


MO2=H2/WO02 — 


COMPUTE MOS 
HS=2*H2=H1 
MOS=HS/WOS 


COMPUTE WEIGHT OF LOG 
V2=R2*ST5+H2*ST4 
VS=RS*ST3+HS*ST6 
WL=V1=R1+VS=RS+2.0* (V2=R2) =WC 
IF (IPRINT-EQe1) GO TO 91 
RETURN 


FOR HeLT»«O = COMPUTE PAYLOAD TYPE=7 
20 CONT INUE 


COMPUTE M02 
ARG1=(R2*STUAST5) #*2=(1604+ST4**2) * ( (R2*ST5) **2=TAK*2) 
IF (ARGieLTe-0.0) GO TO 90 
H2=(=R2*ST5*ST4+SQRT (ARG1) )/(ST4**2+1.0) 
MO2=H2/Ww02 
COMPUTE MO1 
V2=R2*ST5S+tH2*ST4 
DUML= (V2=R2) **2+H2**2 
ARG2=(R1*ST1*(ST2=1.0)) *#2-(1.04+ST1I**2) *( (R1I*(ST2—-1.0) ) **2—-DUM1) 


IF (ARG2eLT.0-0) GO TO 90 
Hi=(=R1*ST1* (ST2—1.0) +SQRT (ARG2))/7(1-0+ST1**2) 
MOL=H1i/WO1 

COMPUTE MOS 
HS=2 6 0*H2=H1 
MOS=HS/WOS 

COMPUTE WEIGHT OF LOG 
VIL=R1I*ST2+HL*ST1L 
VS=RS*ST3+HS*ST6 
WL=V1I"=R1I+VS=RS+20 0* (V2—=R2) =WC 
IF (IPRINT.-EQe1) GO TO 91 
RETURN 


FOR DYeGE.-0 - TAS AT SPAR 1 TYPE=8 


COMPUTE MOS 
ARGI=(RS*ST3*ST6) €#2—(12e0+STEO**2) *( (RS*ST3S) #*¥2=TAS#*2) 
IF (ARG1ieLTe000) GO TO 90 
HS=(=RS*ST3*ST6+SQRT(ARG1) )7(1-¢0+ST6**2) 
MOS=HS/WOS 


COMPUTE MO2 
DUM1=2e0*R1*ST1* (ST2-160) -R2* (ST5=1.0) *ST4-2. 0*HS* (120+ST1¥*2) 
DUM2=(R1I*(ST2-1.0) ) **2+HS* #24 (HS*ST1) *#2=2.0*R1¥*(ST2=1.0) *STI*HS 


al -(R2*(ST5-1.0))**2 


ARG2=DUM1** 2=( 4.0% (ST1**24+1.0)=(ST4**2+1.0) ) *DUM2 

IF (ARG2.LTe0.0) GO TO 90 
H2=(=DUM1+SQRT (ARG2) ) 7(4.0* (ST1¥*2+160)—(ST4*#2+4+1.0) ) 
MO2=H2/w02 


C COMPUTE MO1 
H1=2.0*H2=HS 
MOi=H1/WO1 


C COMPUTE WEIGHT OF LOG 
V1I=R1I*ST2+Hi*ST1 
V2=R2*ST5+H2*ST4 
VS=RS*ST3+HS*ST6 
WL=V1=R1+VS"=RS+260* (V2=R2) -WwC 
IF (CIPRINT«EQe1) GO TO 91 
RETURN 


C FOR DYeLTeO = TAS AT CARRIAGE TYPE=9 


350 CONT INUE 
C COMPUTE MOS 
ARG1=(RS*#ST6*(ST3-1-0) ) *¥2—=(100+ST6**2) * (RSFK2* (STS=1 00) *¥2 
1 =TAS**2) 
IF (ARG1.LT«0-.0) GO TO 90 
HS=(=RS*ST6* (ST3—1.0)+SQRT(ARG1))/(1-0+ST6**2) 
MOS=HS/WOS 
C COMPUTE M02 
DUM1=2.0#R1*ST1*(ST2=1.0) -R2* (ST5-120) *ST4=2. OF HS* (1.04+STL**2) 
DUM2=(R1* (ST2—-1.0) ) **¥24HS**2+( HS*ST1) #*¥2—-2.0%*R1¥* (ST2-1.20) *STI*HS 


1 = (R2*(ST5=1.0).) ** 


ARG2=DUM1**2=(4.0*(ST1*¥*24+1.0) =(ST4**24+1.0) ) *DUM2 
IF (ARG2.LT.0.0) GO TO 90 
H2=(=DUM1+SQRT (ARG2) )/(4.0* (ST1**2+120) -(ST4**2+1.0) ) 
MO2=H2/w02 
C COMPUTE MO1 
H1=2-.0*H2-HS 
MO1=H1/WO1 
C COMPUTE WEIGHT OF LOG 
VI=R1I*ST2+H1I*ST1 
V2=R2*ST5+H2*STY4 
VS=RS*ST3+HS*ST6 
WLEV1=R1+VS=-RSt+2e0* (V2=R2) -WC 
TER INT ee Qs1) GO TO 91 
RETURN 
91 WRITE (OT*#2001) ITYPE* ARG1»ARG29STirST2+ST3SeST4eSTSeST6eLivL2el Se 


1 Ele E2°ESeR1ep R20 RSe DUM1 +s DUM2e H1e H22HSe Vie V29VSeWCe WL 
RETURN 
C 


1 Ele E22 ESeR1eR2eRSe DUM1» DUM2 2 H1 e H2e HSe Vie V20VSeWC 
RETURN 
2000 FORMAT (24H NEG ARG IN FB3 TVPES=PV2nioEtSe7AOClXr Sel S 3% ))) 
2001 FORMAT (6H TYPE=? 12°6E152.7/3(1Xe 8E15-7/) ) 


END 

SUBROUTINE TESTSCITYPE? 1) 
C 
C PROGRAM TO CHECK CONVERGENCE 
(e 


21 


90 WRITE (OT? 2000) ITYPEs+ ARG1PARG2eST1eST2eST3eST4eST5eST6eLirpl2p LS» 


C 


REAL LeLireL2eLS»MO1»MO2* MOS 
COMMON IPRINTs ITER 


COMMON TAr TASeWOL*® WO2 »WOSe Le Se DEFL»# De DY» DSe He DYSe WCe C1 HS1 9 HS2 


COMMON LieL2eLSeE1re E20 ES 


COMMON MO1*MO2»* MOS» XAMO1* X1MO1* XAMO2s X2MO2e XAMOS» XI MOS 


COMMON WL*SG 
INTEGER OT 
OT=61 


- 


WL=ZMOS*WOS*S INH( XAMOS) +MO1*WO1*S INH( XAMO1 ) +26 0*MO2*WO2*S INH( XAMO2 | 


TOL=0.5 

T=2 

IDUM=ITYPE=5 

GO TO (10°20+30*°40) » IDUM 


C PORTION TREATING H.eGE.O0 = ITYPE = 6 


C 


10 


CONT INUE 


C TEST FM7 AND FM8 


FM7=T A=wO1*DY-WO1*MO1* COSH( XAMO1 ) 
FM8=T A=w0O1* DY-wO2*MO02* COSH( XAMO2) 
ie CANS CSO) abitp slat oak 

TE CABS(FM8)=TOE)” 1230139 12 

I=1 


TF .ChiPRis= G10), sGO O50 
TE CIPRINT oEQe1) ‘GO TO) Si 
RETURN 


€ PORTION TREATING HeLT.O = ITYPE 7 


C 
20 


CONT INUE 


FM9=T A=WO1*MO1*COSH( XAMO1 ) =-WO2* (DY=H) 
FM10=TA=WO2 ( DY=H) =wO2*MO2* COSH( XAMO2) 


C TEST FM9 AND FM10 


al 
22 
C 
23 


C 
C 
30 


TE GABSIGEMS)i=1OE)  2ti2ii22 
TES CABSICEMLO) =O) 2Sr2sr22 
T=1 


IF (ITEReEQe10) GO TO 60 
IF (IPRINTeEQe1) GO TO 61 
RETURN 


PROTION TREATING DY.GEe0 - ITYPE 8 


CONT INUE 


C TEST FM11 AND FM12 
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FM11=M01*WO1*COSH(XAMO1) =MO2*WO2* COSH( X AMO2) 
FM12=TAS=wOS*DY=wOS*MOS*COSH( XAMOS) 


te WANSISC aM) SOL) Silo silo 3 


TE .WABSCEML2) —nOl i Sorsomse 
I=1 
IF (ITEReEQe-10) GO TO 70 
IF CIPRINT.EQ.1) GO TO 71 
RETURN 


)- we 


4 
' 


40 


| 
60 
| 61 


70 
ae. 


| C 


C PORTION TREATING DYeLT-O - ITYPE = 9 


CONTINUE 


Se TEST EMIT AND) MLS 


FM11=M01*WO1* COSH( XAMO1) =-MO2*WO2* COSH( XAMO2) 
FM13=T AS@MOS*WOS*COSH( XAMOS) 

IF (ABS(FM11)=TOL) 41°41-42 

IF (ABS(FM13)-TOL) 43°43942 

I=1 


IF (ITEReEQe10) GO TO 80 
IF (IPRINT.EQ.1) GO TO 81 
RETURN 


WRITE (OT?2000) 

CONT INUE 

WRITE (OT*3003) ITER*#FM7eFM8eLirp L2eeLSe Ele E20 ES 

WRITE (OT*#3002) MO1*»MO2»MOSpr XAMO1 » XLMO1 *? XAMO2 2 X2MO2 » XAMOS ep XL MOS 


GO TO 90 

WRITE (OT*2001) 

CONTINUE 

WRITE (OT*#3004) ITER*FMOeFM10°eLirL2eLSeE1eE2rES 

WRITE (OT*#3002) MO1»MO2*+MOSr XAMO1 2 X1MO1»9 XAMO2 X2MO02 0 XAMOS¢ X1MOS 


GO TO 90 

WRITE (OT*#2002) 

CONT INUE 

WRITE (O0T*#3005) ITER*#FM11*eFM12eLireL2reLSreE1sE2eEs 

WRITE (OT*#3002) MO1»MO2+MOS»+ XAMO1+ X1MO1 2» XAMO2’ X2MO20 XAMOS? X1MOS 


GO TO 90 

WRITE (0OT*2003) 

CONTINUE 

WRITE (OT*#35006) ITER® FM11*°FM13eLireLe2reLSe Ele E2rES 

WRITE (0T*#3002) MO1*+MO2»MOSr XAMO1 » X1MO1» XAMO22 X2MO02 » XAMOS* XI1MOS 


CONT INUE 


C COMPUTE CATENARY TENSIONS 


2000 
2001 
2002 
3001 
3002 


3003 
3004 


T1=wO1*MO1*COSH(X1M01) 

T2=wO2*MO02*COSH( X2M02) 

TS=wOS*MOS*COSH(X1MOS) 

V1=wO1*MO1*SINH(XAMO1) 

V2=wO2*MO2*S INH( XAMO2) 

VS=wOS*MOS*S INH( XAMOS) 

Hi=wO01*MO1 

H2=w02*MO2 

HS=wOS*MOS 

WRITE (OT*#3001) Tir T20+>TSe V1ire V2eVSe Hie H2*e HS 

RETURN 

FORMAT (35H TOO MANY ITERATIONS ON HeGE.O CASE) 

FORMAT (35H TOO MANY ITERATIONS ON HeLTeO CASE) 

FORMAT (36H TOO MANY ITERATIONS ON DY«GE.0 CASE) 

FORMAT (29H (T1lrT2*+TSeV1ieV29VSeH1le H2eHS) e 9F11.1) 

FORMAT (50H (MO1*MO2*+MOSe XAMO1? X1MO1 9 XAMO22 X2M02 0 XAMOSP X1MOS) » 
3F8.196F 8.5) 


FORMAT (33H (ITER* FM7e FM8eL1leL2eLSe Ele E2rES) e [52 2F12.6r6F 10601) 
FORMAT (34H (ITER* FM9eFM10eLireL2eLSr Ele E2r ES) ep 15°2F12.6r6F 1061) 


23 


(ee ee eee ee ene ae ee L 


3005 FORMAT (35H (ITER»FM1i1»FM1i2eL1ireL2eLSeE1e E20 ES) oe 152 2F12e696F 1061 


2003 FORMAT (36H TOO MANY ITERATIONS ON DYeLT.O CASE) 
3006 FORMAT (35H (ITER*e FM11°FM13*9L1*sL2eLSeEir E22 ES) ep 150° 2F12.696F100l 


END 


SUBROUTINE DEFL3(ITYPE) 
c é 
C DETERMINATION OF DEFLECTION FOR GIVEN LINE LENGTH AND STATION. 


REAL LeLieL2eLS»MO1»MO2* MOS 
COMMON IPRINT?s ITER 
COMMON TA» TASeWO1L9!WO2%+WOSe Le Se DEFLe De DY e DSe He DYSe WCe C1» HS1 2 HS2 
COMMON LieL2reLSeE1+/E2° ES 
COMMON MO1*MO2*+MOS*e XAMO1 + X1MO1» XAMO2* X2MO02% XAMOS? X1MOS 
COMMON WLe»eSG 
INTEGER OT 
OT=61 
DFINIT=DEFL 
IBO=1 
IBO=0 
IF (IPRINT-EQ.1) IBO=1 
TOL=TA/1000000-0 
i I=0 


C INITIAL ESTIMATE OF DEFLECTION AND PAYLOAD FROM PREVIOUS CASE 
DY=DFINIT*L+S*D 
DYS=DY 
IF (HeGEe0-0) ITYPE=6 
IF (HeLTe0.0) ITYPE=7 
CALL RCONVG(ITYPE) 
WELI=WL 
DF1=DF INIT 
IF (ABS(WL1).eLEe«TOL) RETURN 


C GET SECOND CATENARY SOLUTION 
DFCHG=0.01 
6 DF2=-wWL1/ABS(WL1) *DFCHG+DF1 
IF (DF2.GT.0.0) GO TO 8 
DF CHG=DFCHG*0.5 
GO TO 6 
8 DEFL=DF2 
DY=S*D+L*DEFL 
DYS=DY 
CALL RCONVG( ITYPE) 
WL2=WL 
é IF (ABS(WL2).eLE.TOL) RETURN 
C ENTER SECANT ITERATION LOOP 
10 RAT=1.0 
11 DF3=DF2=( (DF2=DF1) *wL2/(WL2=wL1) )*RAT 
IF (DF3eGTe0.0) GO TO 12 
RAT=RAT/2-0 
GO TO 11 
12 DEFL=DF3 
DY=S*D+L*DEFL 
DYS=DY 
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CALL RCONVG(ITYPE) 
WL3=WL 
IF (TBO«.EQ.1) WRITE (OT*2000) IpWhLieWL2eWL39DF1»DF2°DF3 
IF (ABS(WL3).LEeTOL) RETURN 
I=I+1 
eG iGive 20)! GO. TO 20 
DF1=DF2 
DF2=DF3 
WLi=wL2 
WL2=wL3 
GO TO 10 

20 CONTINUE 
IBO=IBO+i 
IF (IBO~-EQ.2) RETURN 
GO TO 1 

2000 FORMAT (21H ITERATIONS IN DEFL3=» 12? 6F15.6) 

END 


SUBROUTINE CATARM 
C ROUTINE TO GENERATE CATENARY LENGTHS AND MOMENT ARMS 
C 
REAL LeLirL2*>LS»MO1»MO02» MOS 
COMMON IPRINT?* ITER 
COMMON TAs TASe WO1/+WO2* WOSe Le Se DEFL?e De DY 2 DSP He DYSe WCe CL ep HS1 9 HS2 
COMMON LireL2eLSre Ele E29ES 
COMMON MO1*+MO2»*MOSe XAMO1 + X1MO1 » XAMO2 2 X2MO2 9 XAMOS er XI MOS 
COMMON WL?SG 
L1=MO1*FUNC2 (X1MO1» XAMO1) 
L2=M02* FUNC2 ( X2M02 9 XAMO2) 
LS=MOS* FUNC2( X1MOS» XAMOS) 
E1=MO1*FUNC1 (X1MO1»XAMO1) /FUNC2( X1MO1» XAMO1) 
E2=M02*FUNC1 ( X2M02" XAMO2) /FUNC2 ( X2M02» XAMO2) 
~ ES=MOS*FUNC1 (X1MOS»* XAMOS) /FUNC2 ( X1MOS» XAMOS) 
RETURN 
END 


SUBROUTINE CATGEO 
C ROUTINE TO GENERATE CATENARY GEOMETRIC PARAMETERS= PRODUCTION MODEL 


(C 
REAL LeLireL2*eLSeMO1» MO2* MOS 
COMMON IPRINTe ITER 
COMMON TArTASe WOOL) WO2eWOSe Le Se DEFL»? De DY e DSep He DYSo WCe Cie HS1 9» HS2 
COMMON LireLerLSe E11 E2°ES 
COMMON MO1*»MO2*MOS+ XAMO1* X1MO1 9 XAMO20 X2MO2e XAMOS¢ XI MOS 
COMMON WLeSG 
INTEGER OT 
OT=61 

(C 

C ARGUMENT CHECK 
10 IARG=1 

ARG1=D/2-.0/M01 
ARG2=(L=D) /2.0/M02 
ARG3=DS/2-0/MOS 
IF (ABS(ARG1)-GTe87-0) GO TO 90 
IF (ABS(ARG2).GT.87-0) GO TO 90 

C HERE IS WHERE THE PRODUCTION MODEL DIFFERS FROM THE STANDARD 
IF (ABS(ARG3).GTe87.0) MOS=(DS/2.0)/( ARG3/ ABS ( ARG3) *87. 0) 

(e 
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21  X1MO1=FUNC3( DY» MO1+D) 
22 XAMO1=FUNC4 ( DY » M01» D) 
23 X2MO2=FUNC3(DY=He MO2*L=D) 
24  XAMO2=FUNC4(DY—He MO2*L=D) 
C THIS MANUEVER IS TO AVOID UNDERFLOW IN THE CALCULATIONS 
IF (ABS(ARG3).GT.20-0) GO TO 30 
25  X1MOS=FUNC3(DYSr* MOS? DS) 
26  XAMOS=FUNC4(DYS» MOS? DS) 
GO TO 40 
30  X1MOS=DYS/2.0/M0S 
XAMOS==DYS/2.0/MO0S 
C ARGUMENT CHECK 
40 IARG=2 
IF (ABS(XAMO1).G6GT.87.0) GOTO 90 
IF (ABS(X1M01).G6T.87-0) GOTO 90 
IF (ABS(XAMO2).GT.87.0) GOTO 90 
IF (ABS(X2M02) .GT.«8720) GOTO 90 
IF (ABS(XAMOS).GT.87-0) GOTO 930 
IF (ABS(X1MOS).GT.87-0) GOTO 90 
RETURN 
90 WRITE (OT?2000) IARGe ARG1*’ ARG2r ARG3* X1MO1* XAMO1» X2MO22 XAMO2 > 
fl X1MOSr XAMOS 
RETURN 
2000 FORMAT (20H ARG LARGE IN CATGEO?r I4r 9F10.2) 
END 


FUNCTION ASH(X) 

ASH = ALOG(X+SQRT(X*X+1.0) ) 
RETURN 

END 


FUNCTION SINH(X) 

S INH =0.5* (EXP(X) =EXP(=xX) ) 
RETURN 

END 


FUNCTION COSH(X) 

COSH = 0.5*(EXP(X) +EXP(=X) ) 
RETURN 

END 


FUNCTION FUNC3(X+*YrZ) 

FUNCS3 =ASH(XZ(2e0*Y*SINH(2/2.0/Y))) +2/2.0/Y 
RETURN 

END 


FUNCTION FUNC4( XYZ) 

FUNC4 =ASH(X/(2¢0*Y*SINH(Z/2e0/Y))) =-Z2/2.0/Y 
RETURN 

END 


FUNCTION FUNC1(X?ryY) 

FUNC1 = (X*SINH(X) =COSH( X) =Y*S INH(X)+COSH(Y) ) 
RETURN 

END 


FUNCTION FUNC2(XPrY) 
FUNC2 =(SINH(X)=SINH(Y) ) 
RETURN 
END 
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The mission of the PACIFIC NORTHWEST FOREST 
AND RANGE EXPERIMENT STATION is to provide the 
knowledge, technology, and alternatives for present and 
future protection, management, and use of forest, range, and 
related environments. 


Within this overall mission, the Station conducts and 
stimulates research to facilitate and to accelerate progress 
toward the following goals: 


1. Providing safe and efficient technology for inventory, 
protection, and use of resources. 


2. Development and evaluation of alternative methods 
and levels of resource management. 


3. Achievement of optimum sustained resource produc- 
tivity consistent with maintaining a high quality forest 
environment. 


The area of research encompasses Oregon, Washington, 
Alaska, and, in some cases, California, Hawaii, the Western 
States, and the Nation. Results of the research will be made 
available promptly. Project headquarters are at: 


Fairbanks, Alaska Portland, Oregon 
Juneau, Alaska Olympia, Washington 
Bend, Oregon Seattle, Washington 
Corvallis, Oregon Wenatchee, Washington 


La Grande, Oregon 


Mailing address: Pacific Northwest Forest and Range 
Experiment Station 
P.O. Box 3141 
Portland, Oregon 97208 


The FOREST SERVICE.of the U.S. Department of Agriculture 
is dedicated to the principle of multiple use management of the 
Nation’s forest resources for sustained yields of wood, water, 
forage, wildlife, and recreation. Through forestry research, 
cooperation with the States and private forest owners, and 
management of the National Forests and National Grasslands, it 
strives — as directed by Congress — to provide increasingly 
greater service to a growing. Nation. 


